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Abstract. We study scaling relations of compressible isothermal strongly magnetized turbulence 
using numerical simulations with resolution 512^. We find a good correspondence of our results 
with the Fleck (1996) model of compressible hydrodynamic turbulence. In particular, we find 
that the density-weighted velocity, i.e. u = p'/^v, proposed in Kritsuk et al. (2007) obeys the 
Kolmogorov scaling, i.e. Su{k) ~ k^^l^ for the high Mach number turbulence. Similarly, we find 
that the exponents of the third order structure functions for u stay equal to unity for all Mach 
numbers studied. The scaling of higher order correlations obeys the She-Leveque (1994) scalings 
corresponding to the two-dimensional dissipative structures, and this result does not change with 
the Mach number either. In contrast to velocity v which exhibits different scaling parallel and 
perpendicular to the local magnetic field, the scaling of u is similar in both directions. In addition, 
we find that the peaks of density create a hierarchy in which both physical and column densities 
decrease with the scale in accordance to the Fleck (1996) predictions. This hierarchy can be 
related ubiquitous small ionized and neutral structures (SINS) in the interstellar gas. We believe 
that studies of statistics of the column density peaks can provide both consistency check for the 
turbulence velocity studies and insight into supersonic turbulence, when the velocity information is 
not available. 
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INTRODUCTION 

The interstellar medium (ISM) is a highly compressible turbulent, magnetized fluid, 
exhibiting density fluctuations on all observable scales. It has been long realized by 
many researchers that incompressible hydrodynamic, i.e. Kolmogorov, description is 
inadequate for such a medium [see 5, for review] . Scaling relations, if they were obtained 
for the interstellar gas, would be very helpful for addressing many problems, including 
the evolution of molecular clouds and star formation. 

Attempts to include effects of compressibility into the interstellar turbulence descrip- 
tion can be dated as far back as the work by von Weizsacker [20] . There a simple model 
based on a hierarchy of clouds was presented. According to this picture every large cloud 
consists a certain number of smaller clouds, which contain even smaller clouds. For such 
a model von Weizsacker [20] proposed a relation between subsequent levels of hierarchy 

pv/pv-l=(/vAv-l)"^", (1) 

where pv is the average density inside a cloud at level v, Zy is the mean size of that cloud, 
3 is the number of dimensions, and a is constant that reflects the degree of compression 
at each level v. 



The Kolmogorov energy spectrum (~ k^^/'^) follows from the assumption of a con- 
stant specific energy transfer rate e ~ v^/(//v). Lighthill [15] pointed out that, in a com- 
pressible fluid, the volume energy transfer rate is constant in a statistical steady state 

ey = pe-pv2/(//v)=pvV/. (2) 

In an important, but not sufficiently appreciated work. Fleck [6] (henceforth, F96) 

incorporated above hierarchical model with energy transfer in compressible fluid to 
obtain the scaling relations for compressible turbulence. By combining the equations 
(1) and (2) Fleck [6] presented the following set of scaling relations in terms of the 
degree of compression a: 

pi - A^; ^ Ml - VI - (3) 

where Ni and M; are, respectively, the column density of the fluctuation with the scale / 
and the mass of the cloud of size /. The fluctuations of velocities in F96 model entail the 
spectrum of velocities E{k) ~ k~^/^~^^. 

In the spirit of F96 model, [12] proposed to use the density- weighted velocity 
u = p^/^v as a new quantity, for which the Kolmogorov scaling for second order struc- 
ture functions (SFs) can be restored in compressible hydrodynamic turbulence. Their 
hydrodynamic simulations provided the spectrum for u close to -5/3 and they showed 
that in supersonic hydrodynamic turbulence the structure functions of u scale linearly 
with separation. 

Will the F96 model be valid for compressible strongly magnetized XmhuloacQl This is 
the major question that we address in this paper. 



NUMERICAL IVIODELING 

We used an second-order-accurate essentially nonoscillatory (ENO) scheme [see 3, 11, 
for details] to solve the ideal isothermal magnetohydrodynamic (MHD) equations in 
a periodic box with maintaining the V ■ 5 = constraint numerically. We drove the 
turbulence at wave scale k ~ 2.5 (2.5 times smaller than the size of the box) using a 
random solenoidal large-scale driving acceleration. This scale defines the injection scale 
in our models. The rms velocity 5v is maintained to be approximately unity, so that 
V can be viewed as the velocity measured in units of the rms velocity of the system 

and B/ {AnpY^^ as the Alfven velocity in the same units. The time t is in units of the 
large eddy turnover time (~ L/ 5v) and the length in units of L, the scale of the energy 
injection. The magnetic field consists of the uniform background field and a fluctuating 
field: B = 5ext + b. Initially b = and v = 0. We use units in which the Alfven speed 
va =5ext/ (4;rp)^^^ = 1 and p = 1 initially. Structures of density, velocity and magnetic 
field develop completely due to the forcing from uniform initial conditions. The values 
of 5ext have been chosen to be similar to those observed in the ISM turbulence. For 
our calculations we assumed thatfigxt/ (4;rp)^''^ ~ 5B/ {AnpY^'^ ~ 5v. In this case, the 
sound speed is the controlling parameter, and basically two regimes can exist: supersonic 
and subsonic. Note that within our model, supersonic means low j8 (j8 = pgas/ Pmag), 
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FIGURE 1. Spectra of velocity and density-weighted velocity (dashed and sohd lines, respectively) for 
super- and subsonic models (big and small plots, respectively). Spectra are compensated by k^^^. 



i.e. the magnetic pressure dominates, and subsonic means high j8, i.e. the gas pressure 
dominates. 

We present results for selected 3D numerical experiments of compressible MHD 
turbulence with a strong magnetic field for sonic Mach numbers between ^ 0.7 and 
7. The Alfvenic Mach number ~ 0.7. Mach numbers are defined as = 
and = (kl/cA) for the sonic and Alfvenic Mach number, respectively. To study 
effects of magnetization we also performed superAlfvenic experiments with ^a ~ 2. 
All models were calculated with the resolution 512^ up to 6 dynamical times. 

RESULTS 

Kolmogorov Scalings for Supersonic Flows 

In Figure 1 we present the spectra for velocity v and density-weighted velocity u = 
pi/^v for two strongly magnetized models: subsonic (j8 ~ 2) and supersonic (/3 ~ 0.02). 
Naturally, for subsonic model the differences between spectra for v and u are marginal 
and both spectra correspond to Kolmogorov's k~^/^ scaling (see subplot in Fig. 1). 
However, we can see that for the supersonic case, the velocity spectrum gets steeper. 
The steepening corresponds to a ^ 0.23 (from S'v ~ At the same time, the 

spectrum of u matches well the Kolmogorov slope. 

In the original Kolmogorov theory [10, hereafter K41] it was shown that the expo- 

nent of the third order structure function (SF), e.g. structure function of v, Sy (/) = 
(|v(r + /) — v(r) p) ~ l^^, should be equal 1, i.e. ^3 = 1. In Figure 2 we show the SFs 
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FIGURE 2. SFs of the third order for v (dashed line) and u (soHd Hne) compensated by /^^ for 
supersonic MHD turbulence model. Dotted lines correspond to best fitting within the inertial range. Two 
dotted vertical lines bound the intertial range. The sub-panel shows the scaling exponents 1^3 for v (crosses) 
and u (diamonds) as a function of ^^s- 




of the third order for velocity and the density-weighted velocity for supersonic model. 
We checked, that for the subsonic case for both v and u the index ^3 is indeed close 
to unity. For the supersonic case, ^3 increases with ^5 for v, but stays the unity for u 
(see Fig. 2). This suggests that the Kolmogorov universality is preserved for supersonic 
MHD turbulence when density weighting is applied. 



She-Leveque Intermittency Model 

A proper description of turbulence requires higher moments [see 13, for review]. 
Those characterize intermittency, which is in the original K41 model is not ac- 
counted for. A substantial progress in understanding turbulence intermittency is 
related to a discovery by She & Leveque ([19], hereafter SL94), who found a sim- 
ple form for the scaling of exponents i^p of higher order longitudinal correlations 

S^P\l) = (I [v{r + l)-v{r)]-l\P) ~ l^p. While in K41 model = p/3, SL94 provides 
(after modifications by Miiller & Biskamp [17] to more general form) 



C„ = ^(l-x) + (3-D) l-{l-x/{3-D))P/^ 



(4) 



where g is related to the scaling of the velocity v/ ~ /^/^, x is related to the energy 
cascade rate tf^ ~ l^"^ and D is the dimension of the dissipative structures. In hydro- 
dynamic incompressible turbulence, we have g = 3 and x — 2/3. For MHD turbulence 
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FIGURE 3. Scaling exponents for v and u for subsonic (subplot) and supersonic models. The plots 
present unnormalized values of the scaling exponents obtained directly from SFs by fitting the relation 
S^P\l) = al^p within the intertial range, i.e. without using the extended self- similarity [1]. 



the dissipation happens in current sheets, which are two-dimensional dissipative struc- 
tures, corresponding to D = 2 [17]. Thus, for subsonic MHD turbulence we expect 
= I + 1 — (1/3)^/-^ for both velocity and the density-weighted velocity. This is what 
we actually observe in Figure 3 (see subpanel). The same scaling, however, is preserved 
for u for supersonic magnetized turbulence. 



Anisotropics Induced by Magnetic Field 

Magnetic field is known to induce anisotropics of compressible MHD turbulence [see 
9]. Anisotropy increasing with the decrease of scale was predicted for Alfvenic motions 
by Goldreich & Sridhar ([7], henceforth GS95, see also [16]) and confirmed numerically 
for compressible MHD in [3,4]. 

For supersonic motions Figure 4 shows that the scalings for v are much steeper in both 
directions than those predicted by GS95 model (1.215 and 0.874 for || and ± directions 
to the local magnetic field, respectively). However, those slopes still give a close to 
GS95 anisotropy (Zy ~ /|[^'^), which is indicative of the dominance of the Alfvenic 
("incompressible") motions. Note, that in Figure 4 the SFs are obtained in the system of 
reference of the local magnetic field, i.e. the field on the scales of the fluctuations under 

Study. 5 1 ' and S\_ ' denote second order SFs parallel and perpendicular to local magnetic 
field, respectively. 

For u the scalings are significantly shallower (0.882 and 0.744 for || and _L directions 
to the local magnetic field, respectively). The SF in perpendicular direction scales more 
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FIGURE 4. SFs of the second order for v and u in the local reference frame for the supersonic 
experiment. SFs for velocity scale as ~ ;1-215 ^jj^j ^ ^0.874 ^^j. p^j-^Ugj ^jjjj perpendicular directions to 
the local magnetic field, respectively. SFs for u scale as ^ /0 882 ^ ^0.744 ^^j. || j_ directions, 
respectively. 



like incompressible motions, i.e. 5^^^ ~ ll_ . The slope of for u is smaller than the 

corresponding one for v resulting in the reduced degree of anisotropy (/|| ~ /|[-^^^). This 
is well visible for small-scale structures as presented in the Figure 5. Intuitively, this can 
be understood in terms of dense clamps strongly distorting magnetic field as they move 
in respect to magnetized fluid. 



Spectra for Density Fluctuations 

Our results for velocity show that our simulations of strongly magnetized turbulence 
provide a ~ 0.23 for ~ 7. The spectrum of density fluctuations ~ ^^^+^" follows 
from the scaling relation of density (see Eq. 3) according to the F96 model. This suggests 
the existence of a rising spectrum of density fluctuations within the hierarchy of density 
clumps when a > 1/6. 

The Figure 6 shows spectra of the low-value filtered densities for supersonic model. 
Spectra for densities above a threshold of 5% of the maximum value start to be rising. 
The Figure 6 contains also a table of mass and volume filling factors, which shows that 
with increasing threashold, the same amount of matter occupies less space. 




FIGURE 5. Contours of SFs of the second order for v and u (left and right plots, respectively) in the 
local reference frame for the supersonic experiment. The small-scale structures show the reduced degree 
of anisotropy. 




FIGURE 6. Spectra of the low-value filtered densities for strongly magnetized supersonic turbulence for 
three thresholds: 1%, 5% and 10% of the maximum value of p (long dashed, short dashed and dot-dashed 
lines, respectively). Solid line correspond to the spectrum of density without filtering. The embedded table 
shows values of spectral indices and percent of the total mass and total volume occupied by the filtered 
density for subsequent thresholds. 
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FIGURE 7. Scaling relations for the column density (left plot) and for density (right plot) for three 
models of turbulence: subsonic ~ 0.7) and supersonic ^ 2 and 7, see legend). 



Statistics of Column Density Peaks 

We try to make our study more related to observations which usually measure den- 
sities integrated along the line of sight, i.e. column densities, or alternatively study the 
hierarchy of observed clump masses (see Eq. 3). F96 model assumes the existence of an 
infinitely extended hierarchy. In our computations the structures are generated by turbu- 
lence at scales less than the scale of the computational box. Therefore the F96 scaling 
relations (in Eq. 3) should be modified as follows 

Our procedure of obtaining the scaling relation from column density maps is similar 
to that in [12], with the difference that they dealt with 3D data, while we deal with 
2D data. First, we seek for a local maximum of column density. Then we calculate the 
average column density within concentric boxes with gradually increasing the size /. In 
case of determining the relation for A//, instead of averaging we apply the integration 
over the boxes. Naturally, the results should correspond to each other. 

In Figure 7 we present an example of scaling relation for column density and density 
for three models of turbulence with ~ 0.7, 2, and 7. One can note, that the relation 
becomes more steep with the sonic Mach number within the intertial range. The fractal 
dimensions can be calculated from the relation D„, = 3 + 7 [see 12], where 7 = — 3a 
is a slope estimated from the plot within the intertial range. For our models the fractal 
dimension ranges from D,„ ~ 2.5 for the highly supersonic models to ~ 2.9 for 
subsonic model. Respectively, the compressibility coefficients for presented models 
a ^ 0.04 for ~ 0.7 to a ~ 0.19 for ~ 7. The latter roughly consistent with 
a obtained for the velocity SF measurement in Section "Kolmogorov Scalings for 
Supersonic Flows". The differences are probably due to insufficient statistics of rather 
rare high density peaks. In general, the filling factor of a peak decreases with the 
maximum density of this peak, which means that the higher maximum of density the 
peak has the smaller space it occupies. 



Variations of Scalings Induced by Fluid Magnetization 



What is the effect of magnetic field on the M-scaling? The spectra, third and higher mo- 
ments of correlations obtained for our superAlfvenic simulations with ~ 2 happen to 
be very similar to the case of strongly magnetized turbulence. Our results indicate that, 
unlike velocity, u is much less affected by magnetic field. Naturally, in super-Alfvenic 
turbulence the anisotropies induced by magnetic field are not observed at larger scales 
within the intertial range (cf. the last paragraph of Section "Astrophysical Implications". 

ASTROPHYSICAL IMPLICATIONS 

Dependence of a on the extend of inertial range. If we combine several facts 
together, namely, (a) that a is a function of Mach number, (b) that the maximum of 
density correspond to the dissipation scale, e.g. shock thickness scale Ijisg, (c) that 
the amplitude of density in peaks scales as the mean density times we have to 
conclude that as the inertial range from the injection scale /,„y to l^iss increases, for a 
given Mach number, a should decrease. Connecting these facts we get the following 
relations, ppeak ~ ~ {Unj/ldiss)^"', which gives the dependence of a on knj/ldiss 
and namely, a ~ Xog^Ms / ^og{linj / Idiss) ■ An interesting consequence of this would 
be a prediction of Kolmogorov scaling for supersonic velocities when the injection and 
dissipation ranges are infinitely separated. Consequently, the steeper velocity spectra 
reported in [18] can be interpreted as an indication of a limited inertial range. Further 
research justifying such conclusions is required, however. 

SINS of supersonic turbulence. Ubiquitous small ionized and neutral structures 
(SINS) are observed in interstellar medium [see 8]. Their nature is extremely puzzling 
if one thinks in terms of Kolmogorov scalings for density fluctuations. The fact that the 
spectrum of fluctuations of density in supersonic turbulence is shallower that the Kol- 
mogorov one is well-known [see 11, and references therein]. However, just the differ- 
ence in slope cannot explain the really dramatic variations in column densities observed. 
The present paper provides a different outlook at the problem of SINS. We see that, 
while low amplitude density fluctuations exhibit Kolmogorov scaling [2, 11], high peaks 
of density correspond to a rising spectrum of fluctuations. Thus, observing supersonic 
turbulence at small scales, we shall most frequently observe small amplitude fluctua- 
tions corresponding to Kolmogorov-like spectrum of density fluctuations. Occasionally, 
but inevitably, one will encounter isolated high density peaks. An alternative mechanism 
for getting infrequent large density fluctuations over small scales is presented in [14] and 
is related to current sheets in viscosity-damped regime of MHD turbulence. 

Clumps and star formation. Interstellar medium is known to be clumpy. Frequently 
clumps in molecular clouds are associated with the action of gravity. Our study shows 
that supersonic turbulence tend to produce small very dense clumps. If such clumps 
happen to get Jean's mass, they can form stars. Therefore, star formation is inevitable in 
supersonic turbulence. However, the efficiency of star formation is expected to be low, 
as the filling factor of peaks decreases with the increase of the peak height. Inhibiting 



of star formation via shearing may dominate in terms of influencing of star-formation 
efficiencies. 

SUMMARY 

In paper above we have studied the scaUng of supersonic MHD turbulence. We found 
that: 

• Fleck 1996 model is applicable to strongly magnetized compressible turbulence. 

• Spectra and structure functions of density-weighted velocities are consistent with 
predictions of the Kolmogorov theory. 

• Intermittency of density-weighted velocity can be well described by the She- 
Leveque model with the dimension of dissipative structures equal 2 [17]. 

• Strongly magnetized supersonic turbulence demonstrate lower degree of anisotropy 
if described using the density-weighted velocity. 

• The high peaks of column densities exhibit increase of the mean values of column 
densities with the decrease of scale, which may be relevant to the explanation of 
SINS. 
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